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1 Introduction 

Rare B meson decays, induced by the flavor-changing neutral current (FCNC) b ^ s (d) 
transitions is one of the most promising research area in particle physics. Theoretical 
interest to the rare B decays lies in their role as a potential precision testing ground for 
the Standard Model (SM) at loop level. Experimentally, these decays will provide a more 
precise determination of the elements of the Cabibbo-Kobayashi-Maskawa matrix (CKM), 
such as Vtq, (g = d, s, b) and Vub and CP violation. 

The impressive experimental search for the study of the B meson decay will be carried 
out in future, when new experimental facilities, especially the i?-factories at Belle ||l| and 
BaBar [^, are upgraded, and after which the large number (10®-10^^) of B hadrons that 
is expected to be produced in these factories, will allow measuring the FCNC decays of B 
mesons. 

In the first hand, the most reliable quantitative test of FCNC in B meson decays is 
expected to be measured in the B -^ Xs(d)^~^^~ decay. The matrix elements of the b — *• 
si~^£~ transition contains terms describing the virtual effects by ti, cc and uu loops which 
are proportional to combination of the CKM elements Vi^V^*, VcbV*^ and VubV*^ respectively. 
Using the unitarity condition of the CKM matrix and neglecting Vu^Vus i^ comparison to 
VtbV^l and VcbV*^, it is obvious that the matrix element for the b -^ si~^i~ decay involves 
only one independent CKM matrix factor, VtbV^*, so that CP-violation in this channel is 
strongly suppressed in the SM. 

The situation is totally different for the b — ^ di~^i~ decay, since all three CKM factors 
are of the same order in SM, and therefore can induce considerable CP violation in the 
decay rate difference of the b -^ dC.^C.~ and b -^ di^i~ processes (for the current status of 
B -^ Xs(rf)£"'"£~ decay in SM, see 0] and the references therein). So, the b -^ di~^i~ is a 
promising decay for establishing CP violation in B mesons. 

The rare B meson decays are also very sensitive to the 'new physics' beyond SM, such 
as the two Higgs doublet model (2HDM), minimal supersymmetric extension of the SM 
(MSSM) 0, etc. 

One of the most popular extension of the SM is the 2HDM 0, which contains two 
complex Higgs doublets rather than one, as is the case in the SM. In the 2HDM, the FCNC 
that appear at the tree level, are avoided by imposing an ad hoc discrete symmetry 0. 
One possible way to avoid these unwanted FCNC at tree level, is to couple all fermions to 
only one of the two Higgs doublets (Model I). The other possibility is the coupling of the 
up and down quarks to the first and second doublets, respectively (Model II). 

Models I and II have been extensively investigated theoretically and tested experimen- 
tally (see [^] and references therein). The 2HDM without the ad hoc discrete symmetry 
was analyzed in 0-^1 . It is clear that the tree level FCNC appears in this model, however 
their couplings involving first and second generations must be strongly suppressed. This 
conclusion is the result of the analysis of low energy experiments. Therefore, Model HI 
should be parametrized in a way that suppresses the tree level FCNC couplings of the first 
generation, while the tree level FCNC couplings involving the third generation can be made 
non-zero as long as they do not violate the existing experimental data, i.e., B-B^ mixing. 

In this work, following 0, we assume that all tree level FCNC couplings are negligible. 
However, even under this assumption, the couplings of fermions with Higgs bosons may 



have a complex phase e (see |0 and |T^). The constraints on the phase angle 6 in the 
product XttKb (see below) of Higgs-fermion couplings imposed by the neutron electric dipole 
moment, B-B^ mixing, po parameter and b ^ sj decay are discussed in [|^. 

The aim of the present work is the quantitative investigation of the CP violation in 
the inclusive B -^ X^ i~^i~ decay in context of the general 2HDM, in which a new phase 
parameter is present (see below). In other words, this model contains a new source of CP 
violation whose interference with the SM phase can induce considerable difference in the 
CP violation predicted by the SM. To find an answer to the question of, "to what extend 
the new physics effects the results of the SM" , is the main goal of the present work. 

The paper is organized as follows. In Sect. 2, we present the necessary theoretical 
background for the general 2HDM and calculate the branching ratio, CP violation and 
forward-backward asymmetry in the b — ^ di~^i~ decay. Finally, Sect. 3 is devoted to the 
numerical analysis and concluding remarks. 

2 The formalism 

Before presenting the necessary theoretical expressions for studying b -^ di~^i~ decay, let 
us briefly remind the main essential points of the Model III. In this model, without loss of 
generality, we can choose a basis such that the first Higgs doublet creates all fermion and 
gauge boson masses, whose vacuum expectation values are 



(^i) 



/ \ 



V x/2 y 
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In this basis the first doublet Hi is the same as in the SM, and all new Higgs bosons result 
from the second doublet H2, which can be written as 

1 / V2G+ \ ^ I V2H+ 

where G~^ and G° are the Goldstone bosons. The neutral Xi ^^^ X2 ^^^ ^^^ ^^^ physical 
basis, but their linear combination gives the physical neutral H^ and h^ Higgs bosons: 

Xi = H^ cos a — h^ sin a , 
X2 = H^ sin a + h^ cos a . 

The general Yukawa Lagrangian can be written as 

Cy = vI'AlHiU.r + rif^Q^LH{DjR + ^f^Q^LH2UjR + ^f^Q,LH2V,R + h.c. , (1) 

where i, j are the generation indices, ifi,2 = '^cr2-f^i,2, Vij ^^^ iij ; i^ general, are the 
non-diagonal coupling matrices, Qn is the left-handed fermion doublet, Ur and Dr are 
the right-handed singlets. In Eq. (1) all states are weak states, that can be transformed to 



the mass eigenstates by rotation. After this rotation is performed, the Yukawa Lagrangian 
takes the following form (only the part of the Yukawa Lagrangian that is relevant to our 
analysis) 



C-Y = -H+U [VcKueR - rVcKML\ V . (2) 

The FCNC couplings are contained in the matrices S,^ '■^. In the present analysis, we will 
use a simple ansatz for ^^ ''^ |H|, 






y2m 



w 



This ansatz guarantees that FCNC for the first generations are strongly suppressed since 
it is proportional to the small quark mass. It follows from Eq. (3) that, we can safely 
neglect the neutral Higgs boson exchange diagrams, and only charged Higgs boson gives 
new contribution to the h -^ di^i~ decay. Note that Xij are complex parameters of order 
C(l) (see 0), i.e., A^- = |Ajj|e*^. Being complex, Xij allow the charged Higgs boson to 
interfere destructively or constructively to the SM results. In other words, branching ratio, 
as well as CP asymmetry, can get increased or decreased for the b -^ di~^i~ decay in Model 
III. For simplicity we choose ^'^'■^ to be diagonal to suppress all tree level FCNC couplings, 
and as a result Xij are also diagonal but remain complex. Note that the results for Model 
I and Model II can be obtained from Model III by the following substitutions: 

Xtt = cot P Xbb = — cot P for Model I , 

Xu = cot(3 Afcb = + tan /3 for Model II . (4) 

After these preliminary remarks, let us return our attention to the b — > d£~^i~ decay. 
The powerful framework into which the perturbative QCD corrections to the physical decay 
amplitude incorporated in a systematic way, is the effective Hamiltonian method. In this 
approach, the heavy degrees of freedom in the present case, i.e., t quark, W^, H^, h^, H^ 
are all integrated out. The procedure is to match the full theory with the effective theory 
at high scale /i = mw, and then calculate the Wilson coefficients at lower /i ~ C(mb) 
using the renormalization group equations. In our calculations we choose the higher scale 
as /i = mw, since In the version of the 2IIDM we consider in this work, the charged Higgs 
boson the charged Higgs boson is heavy enough {mH± > 210 GeV see [|1T|) to neglect the 
evolution from mH± to rriw. In the version of the 2HDM we consider in this work, the 
charged Higgs boson exchange diagrams do not produce new operators and the operator 
set is the same as the one used for the b — > di'^£~ decay in the SM, but the values of the 
Wilson coefficients are changed at mw scale. The effective Hamiltonian for the b — *■ di~^i~ 
decay is The effective Hamiltonian for the b — > di^i~ decay is 



12, 13 
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where 






and Ci are the Wilson coefficients. The exphcit form of all operators Oj can be found in 



T^, |13] . Using the effective Hamiltonian, the matrix element of the h — > M'^f. decay takes 



the following form: 



Gpa ^^ ^^^1 ^e//^ 



- 2Cf^'^dta^,p''{l+^5)biY^\ , (5) 

where g^ = [pi + ^2)^ is the invariant dilepton mass squared, pi and p2 are the four- 
momentum of leptons. In Eq. (5) all Wilson coefficients are evaluated at the n = mi, scale. 
As has already been noted earlier, in the model under consideration the charged Higgs 
boson contributions to leading order at /x = mw scale modify only values of the Wilson 
coefficients, i.e.. 
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The coefficients Cf^^^'\m,w) to the leading order are given by 
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and sin^ 6'vK = 0.23 is the Weinberg angle. The coefficient Cj^dj,) at the scale /i = 0(mb) 
in next to leading order, taking into account the charged Higgs boson contributions, is 



calculated in 11 



47r 

where Cj{mb) and Cj''^ (mb) are the leading and next to leading order contributions, whose 
explicit forms can be found in [|T^]. In our case, the expressions for these coefficients can 
be obtained from the results of [|ll|] by making the following replacements: 

|r|'^|A«|' and XY*^\XuXtb\e'' . 

In the SM, the QCD corrected Wilson coefficient Cglrriiy), which enters to the decay am- 



plitude up to the next leading order has been calculated in [12, 13 1. In order to calculate 

in |[T^. Hence, 



2HDM 



at /i = mfe scale, it is enough to replace Cg (mw) by C, 



■f2HDM 
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including the next to leading order QCD corrections, Cg{mb) can be written as: 



cM = cr^'\^^) 



1 H oj(s) 
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3Ci(/i) + C2ifi)] - ^giO, s) [Csifi] + 3C,ifi) 



^g{l, s) [4C3(^) + 4C4(/i) + 3CM + CM] - ^g{0, s) [CM + 3C,{^) 
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3Csifx) + C,ifx) + 3C,{fx) + Ce{fx) 



(10) 



where m^ = mjmb , s = p'^/ml C(°)(/i) = 3Ci{fi)+C2{fi)+3C3{fi)+C^{fi)+3C5{fi)+CQ{fi) 
and 
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represents the O (a^) correction from the one gluon exchange in the matrix element of Og, 
while the function g {rhc, s) arises from one loop contributions of the four-quark operators 



Oi-Oq, whose form is 






+ J\l-yMQil-y, 
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where j/j = Amf/p^. The Wilson coefficient Cio does not receive any new contribution in 



evolution from /i = rriw to fi = rrib scale, i.e., Cio(wib) = C 



2HDM, 
10 ' 



mw)- The Wilson coeffi- 



cients Cg receives also long distance contributions, which have their origin in the real mm, dd 
and cc intermediate states, i.e., p, u and J/-?/' family. These contributions must be added 
to the complete perturbative results. In the current literature, there exist four different 
approaches in taking into account cc resonance contributions: a) HQET based approach 
|T^, b) the AMM approach |T^, c) the LSW approach [|I^, and d) KS approach [|1^]. In 
the present article, we choose the AMM approach, in which these resonance contributions 
are parametrized using a Breit-Wigner shape with the normalization fixed by data. The 
effective coefficient Cg including the p, u and J /if) resonances are defined as 



a 



eff 



Cg(/i) + Yresis) 



where Yres in NDR scheme is given by 



Yr. ^ -|4(--(.)^A.[3.,p)..,.)]) E,i^^^^ 



- A.,(m., .) [3C,(p) + CM\ ^g^ M^^ - g^ - .M^T J " 

Moreover, the experimental data determines only the product kC^'^^ = 0.875 [l^, which is 
kept fixed. 

Using Eq. (5), the double differential decay rate can be calculated straightforwardly. 
Neglecting the lepton masses and performing summation over final leptons and d quark 
polarizations, the double differential branching ratio takes the following form (the masses 
of the leptons and d quark are all neglected) : 
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- 4^4(5, z) [Re{C'j"iC'g^^ - CioY) + Re{C',"iC'g" + Cio 
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where s = q^ /ml and z = cos 9, and 
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In deriving the above expression, we have normahzed the branching ratio to the branching 
ratio of the semileptonic b -^ civ decay which has been measured experimentally, in order 
to be free of the uncertainties coming from h quark mass. The normalization factor Bq is 
given as 



So = 5 



SL 



3«' \VtdV:,\^ 
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and the phase factor fifhc), and the 0{as) QCD corrected factor |T^ K,{mc) of the b — >• civ 
decay are given by. 
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After integrating over z in Eq. (13), we get 
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(15) 



In order to avoid uncertainties arising from long distance effects, we shall work above the 
p, uj and below the J /if) resonance regions, in the so-called low-g^ region of 

1 GeV'^ <q^ <Q GeV^ , 

and all further numerical analysis will be performed in this region of q^. In principle, 
higher resonance states like p', u', are all expected to contribute to the total branching 
ratio. However their branching ratios are relatively small, and hence are neglected. Per- 
forming integration over g^ in the above-mentioned region, we obtain the partly integrated 
branching ratio 



AB 



6/m'^ 



b ,.dB(b 
ds 
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together with AB for the CP-conjugate decays b -^ dl'^l , and the branching ratio averaged 
over the charge-conjugated states 
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The CP asymmetry for the b 
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and further substituting A^ -^ A* for the conjugated process b -^ di~^i , for the CP 
asymmetry in the partial rate we get, 
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We also have studied the integrated forward-backward asymmetry, whose definition is as 
follows 
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FB 
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3 Numerical analysis 



The values of the main input parameters we have used in the numerical analysis are as 
follows: sin^Ow = 0.2255, m* = 173.8 GeV, a = 1/129, Bsl = 0.104. Further, the 
Wolfenstein parametrization of the CKM matrix |18[ with A = 0.819 and A = 0.2196 [1^ has 
been used. In this parametrization Vcb = AX^, VtdVth = A\^{l—p+iff), V*^Vud = AX^{p—iff), 
where p = p(l — A^/2) and f/ = ri{l — A^/2). Fits of the CKM matrix elements were 
performed in |2^. In further analysis p = 0.3, rj = 0.34 have been used. For the values 
of the parameters \Xtt\ and \Xbb\, we have used the results given in |0], i.e., |Ati| < 0.3 and 
\Xbb\ = 50. 

In Fig. 1 we present the dependence of the partly integrated (i.e., in the region 1 GeV^ < 
q"^ < Q GeV"^) branching ratio for the b — >■ d£~^£~ decay on the phase angle 6 and the charged 
Higgs boson mass mH±, in units of iV^d^fc/VtsV^^I . It is clearly seen that the value of the 
branching ratio varies in the range 2.1 x 10~^ -7- 2.3 x 10~^, and attains at its maximum 
value at ^ = TT. The same figure also depicts that, the contributions from SM and the 
charged Higgs boson interfere constructively in the region < ^ < tt, while destructively 
in the region it < 6 < 2tt. The surface graphics in Fig. 2 illustrates the dependence of the 
partly integrated forward-backward asymmetry Aps on the phase angle 9 and the charged 
Higgs boson mass mH± ■ The forward-backward asymmetry is independent of the value of 
the charged Higgs boson at points 9 = 7r/2 and 9 = 37r/2, where Aps = 0. The value of 



8 



the forward-backward asymmetry ranges between the values — 0.1 -i- +0.1 and its maximum 
value is at ^ = vr. In Fig. 3 we give the dependence of the partly integrated CP violating 
asymmetry on the phase angle 9 and the charged Higgs boson mass mH± ■ We observe that 
the CP violating asymmetry varies in the range —0.038 -. — 0.058 and gets its maximum 
value 6 = TT. As the mass of the charged Higgs increases the range of variation of Acp 
decreases. For completeness we present the SM prediction for the branching ratio and Aqp 
values too. In AMM approach, SM predicts the following values for the branching ratio 
and Acp: AB = 9.61 x 10~^, Acp = 0.044 ^ (we have used the same values for the input 
parameters to be able to compare the predictions of SM and Model III). If we compare our 
results of AB and Acp with that of the SM predictions, we observe that AB in both models 
are very close to each other, but the value of Acp is quite different, keeping in mind that, 
in the present work the central values of all input parameters have been used. Therefore 
our conclusion is that, in establishing Model III along the lines considered in this article, 
Aqp is more efficient than AB. 
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Figure captions 

Fig. 1 The dependence of the partly integrated branching ratio of the b -^ di~^i~ decay on 
the mass of the charged Higgs boson rriHi and the phase angle 9, in units of l^ttiV'^b/VtsV^^I . 

Fig. 2 The dependence of the partly integrated forward-backward asymmetry on the 
mass of the charged Higgs boson mH± and the phase angle 9. 

Fig. 3 The same as in Fig. 2, but for the CP violating asymmetry. 
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